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ABSTRACT 

A method for determining the variations of the resistivity index (Rl) of a family of 
rock samples of complex pore structure with at least two pore networks, as a 
function of the water saturation (Sw), and in the presence of a non-oil conducting 
fluid. The volume fraction (f 1f h, h) occupied by each pore network and the 
distribution of the pore throats in the various pore networks are determined for 
each sample of the family. The values of coefficients relating the variation of its 
electrical resistivity as a function of its water saturation (Sw) are experimentally 
determined on a sample of the family used as a reference sample. The resistivity 
index (Rl) of all the samples of the family is then determined on the basis of the 
variation of parameters describing the layout of the pore network, and using the 
values of the coefficients measured on the reference sample. The method allows 
measurement of a continuous resistivity index curve which is not obtained, as in 
conventional techniques, from only a limited number of points at capillary 
equilibrium. 



BACKGROUND OF THE INVENTION 

Field of the Invention 

[0001] The present invention relates to a method for determining the 
resistivity index Rl, as a function of the water saturation, of certain rocks of 
complex porosity. 

Description of the Prior Art 

[0002] Evaluation of carbonate reservoirs is a particularly difficult task for 
petrophysicists, who still lack precise knowledge concerning the carrying 
properties within these porous media. In relation to siliclastic rocks, carbonates 
may be simpler on the mineralogical plane, but they are incomparably more 
complex in terms of pore structure and surface properties. The largely biological 
origin of the sediments, combined with various diagenetic processes, leads to 
complex pore structures that may be very different from one reservoir to the 
next. For many carbonate systems, resistivity data calibrations carried out in the 
laboratory contradict the observations made in the field (anhydrous production, 
capillary pressure) and the direct water saturation measurements performed on 
preserved cores. 



[0003] The prior art is notably defined by the following publications : 

- Bouvier L. et al., Reconciliation of Log and Laboratory Derived Irreducible 
Water Saturation in a Double Porosity Reservoir, Advances in Core Evaluation, 
edited by Worthington and Longeron, Gordon and Breach Science Publishers, 

- Dixon J.R. et al. (1990), The Effect of Bimodal Pore Size Distribution on 
Electrical Properties of some Middle Eastern Limestone, Soc. Petr. Eng. 20601, 
7 th SPE Middle Eastern Oil Show, Bahrain, pp.743-750, 

- Fleury M. (1998), TRIM : a Fast Resistivity Index Measurement Method", 
Proceedings of the International Symposium of the Society of the Core 
Analysts, Den Hague, 

- Fleury M. et al. (2000) "Frequency Effect on Resistivity Index Curves Using 
a New Method", Proceedings of the 41 st Annual SPWLA Symposium, Dallas, 

- Moore C.H. (2001) "Carbonate Reservoirs, Porosity Evaluation and 
Diagenesis in a Sequence Stratigraphic Framework", Developments in 
Sedimentology 55, Elsevier Editions, 

- Petricola M.J.C. et al. (1995) "Effect of Microporosity in Carbonates: 
Introduction of a Versatile Saturation Equation", Soc. Petr. Eng. 29841, SPE 
Middle Eastern Oil Show, Bahrain, pp.607-615, 

- Sen P.N. et al. (1997) "Resistivity of Partially Saturated Rocks with 
Microporosity", Geophysics, Vol.62, No.2, pp.41 5-425. 
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[0004] Understanding and prediction of the effect of the structure of the pore 
network, of the wettability and of the electrical properties of carbone rocks is a 
real scientific challenge, theoretically as well as experimentally. In fact, correct 
evaluation of these parameters has a major impact on estimation of the oil in 
place, notably for Middle Eastern giant fields, because the difference in relation 
to the standard values of Archie's exponents m and n in relation to value 2 is so 
great that estimation of the water saturation can vary by more than 20 %. 

[0005] Various experimental observations have shown the existence of 
distinct pore populations : micropores, macropores and mesopores, with 
different coexistence degrees. In general, the resistivity index curve 
RI(Sw)=Rt(Sw)/Ro, where Rt is the resistivity of the rock to a water saturation 
Sw, and Ro the resistivity for Sw=1, cannot be described by a power law 
(second Archie's law RI=Sw~ n ), that is n is a function of the saturation itself. 
Microporosity can act as a parallel path for the current, which leads to a 
decrease in the values of n and, therefore, to a gradual insensitivity of the 
resistivity to saturation, as observed on clayey sands. Microporosity can also be 
the cause of the low measured values of n (typically 1 .45). It has also been 
observed that n can increase considerably under certain conditions, and there 
might be a connection between curve RI(Sw) and the capillary pressure curve. 
The increase in the values of n is also a known effect of the wettability, which 
tends to favor aqueous phase discontinuity and therefore to increase the 
resistivity to water wettability. The effect of the wettability can lead to either a 
sudden increase of n, or to a high value n without discontinuity, which can lead 
to confusion. 
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[0006] The experimental curves RI(Sw) which have been drawn from various 
prior works, already mentioned above, and from observations can have (Fig.1) 
four distinct shapes which do not always meet Archie's laws : 

- type I : can be typical of carbonates from the Thamama formation, 

- type II : straightens at intermediate saturation and flattens at low saturation 
(present study), 

- type III : single slope at low saturation, extrapolation at Sw=1 above lr=1 , 
and 

- type IV : typical of oil-wet systems, high values of n that can increase further 
at low saturation. This is also valid for clasts. 

[0007] It can thus be seen that, in the presence of rocks of complex porosity 
which do not meet Archie's laws, a large number of costly measurements is 
necessary to take account of the variability of the pore structure. 

SUMMARY OF THE INVENTION 

[0008] The method according to the invention determines the variations of 
the resistivity index (Rl) of a family of rock samples of complex pore structure as 
a function of the water saturation (Sw), in the presence of a non-conducting 
fluid. The method of the invention comprises the following stages : 

- for each sample of the family comprising at least a first and a second pore 
network, determining a volume fraction occupied by each pore network by 
applying to the various samples an NMR type relaxometry technique; 
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- for each sample of the family, measuring by mercury injection the pore 
throat distribution in the various pore networks; 

- determining experimentally on a sample at least of the family used as the 
reference the values of coefficients relating the variation of the sample's 
electrical resistivity as a function of the sample's water saturation; and 

- determining the resistivity index (Rl) of all the samples of the family on the 
basis of the variation of parameters describing the layout of the pore network 
and using the values of the coefficients measured on the reference sample. 

[0009] According to a first implementation mode suited for samples 
comprising two pore networks, the values of the coefficients relating a total 
conductivity of the sample to a conductivity of the first and second pore 
networks and to respective water saturations of the two networks are 
determined from the reference sample. The resistivity index is calculated from 
the respective volume fractions of the two pore networks and from the value of 
the mean saturation from which the network with the smaller pores is invaded 
by the non-conducting fluid. 

[0010] According to a second implementation mode suited for samples 
comprising a third pore network, the values of coefficients relating a total 
conductivity of a sample to the conductivity of the first and second pore 
networks and to the respective water saturations of the first two networks are 
determined from the reference sample. The resistivity index is calculated from 
the respective volume fractions of the three pore networks, from a value of a 
mean saturation from which the network having the smaller pores among the 
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first two pore networks is invaded by the non-conducting fluid, and from the 
value of the mean saturation from which the network having the larger pores 
among the first two pore networks is invaded by the non-conducting fluid. 

[0011] The resistivity index variations obtained by applying the method are in 
best agreement with the experimental curves whose examples are illustrated in 
Fig.1. The method is based on a powerful experimental technique allowing 
measurement of a continuous resistivity index curve, unlike the conventional 
technique which describes curve Rl only by a limited number of points at 
capillary equilibrium. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0012] Other features and advantages of the method and of the device 
according to the invention will be clear from reading the description hereafter, 
with reference to the accompanying drawings wherein : 

- Fig. 1 shows a diagrammatic representation, in logarithmic coordinates, of 
various types of curves Rl as a function of the water saturation Sw observed 
experimentally with curves I to III being due to pore structure and curve IV is 
typical of the effects of wettability which is not specific to carbonates; 

- Fig. 2 shows a pore size distribution of three carbonates studied : (reservoir 
carbonate, Estaillade limestone and Brauvillier limestone) obtained by NMR 
with the pore volume fraction being to the left of the vertical segment; 

- Figs. 3a and 3b show the sensitivity of the model with two pore networks 
(DPC) to the saturation Sc at which the micropores are invaded by oil (3a) and 
to a which is the initial conductivity ratio of the two pore populations (3b); 
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- Figs. 4a and 4b show the application of model DPC to sample RC at 
ambient temperature (4a) and under the reservoir conditions with dead oil 
(4b) with the experimental curve being drawn as a thin line and the model being 
drawn as a thick line with parameters n-i, n 2 and a being adjusted, fi and Sc 
being measured; 

- Fig. 5 shows an invasion scenario in the model with three pore networks 
(TPC); 

- Fig. 6 shows an electrical layout modelling of the three networks in the TPC 
model; 

- Fig. 7 shows a measured (thin line) and model (thick line) curve Rl for 
sample EL (Estaillade limestone); 

- Fig. 8 shows a measured (thin line) and model (thick line) curve Rl for 
sample BL (Brauvillier limestone); and 

- Fig. 9 shows a simplified calibration sequence for the resistivity logs with for 
non-Archie type rocks, Sw(RI) is necessary and the laboratory data generally 
providing relation RI(Sw). 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

[0013] A number of samples of the same family exist for which calculation of 
the resistivity index of each sample is desired. The method described hereafter 
allows these indices to be obtained, notably from : 
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- experimental measurements performed on each sample of the family, giving 
the number of pore networks of different size and the volume fraction occupied 
by each pore network, as well as the pore throat distribution for each one, and 

- an experimental measurement performed on a single sample of the family, 
used as a reference sample, of the variation of its electrical resistivity as a 
function of its water saturation (Sw). 

[0014] Nomenclature 

Cto = total conductivity at Sw=1 

Q = conductivity of the i-th population at Swi=1 

fi = NMR volume fraction of pores of the i-th population, i=1 ,2,3 

Rl = resistivity index, Rt/RoSw = mean water saturation 

Swj = water saturation of the i-th pore population 

Sc = mean saturation at which the micropores are invaded 

Sm = mean saturation at which the macropores are invaded 

cci i2 = conductivity ratio C1/C2 and C1/C3 respectively. 

[0015] Experimental data on samples 

NMR relaxation time measurements 

[0016] The various pore networks within the rock samples and the volume 
fraction occupied by each network are first determined by means of a well- 
known NMR type relaxometry method. An example of such an NMR 
relaxometry method is for example described in French patent application EN- 
02/... filed by the applicant. 
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[0017] We present the measurements performed on three selected 
carbonate samples. By applying an NMR type known relaxometry method, it is 
seen that all the samples are characterized by a double or triple pore size 
distribution. The first one is a reservoir carbonate sample from the Middle-East, 
classified as a packstone. It was cleaned by means of different solvents at high 
temperature before being subjected to the measurements. Under ambient 
conditions, the experiment was carried out with refined oil and synthetic 
reservoir brine. The experiment carried out under reservoir conditions was 
conducted with filtered degassed reservoir oil at the reservoir pressure. After 
cleaning, the wettability of the sample is characterized by a moderate water 
wettability and, after aging, by a high oil wettability. The other two samples are 
outcrop carbonates (water-wet). The porosity of the Brauvillier limestone (BL) is 
essentially intergranular and due to the oolite cortex. The porosity of the 
Estaillade limestone is both inter and intragranular. The pore populations are 
separated by a factor 10 at least. 

Resistivity index measurements 

[0018] The resistivity index of the samples is then measured. The fast 
resistivity index measurement method (FRIM) described in patent EP-974,839 
filed by the assignee can be advantageously used. 

[0019] According to this method, a forced oil-water displacement is carried 
out on a small core of length 2.5 cm and diameter 4 cm for example. This 
displacement is very close to a porous plate displacement method, except that 
capillary equilibrium is not necessary. Under ambient conditions, only two or 
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three pressure stages are used. Under reservoir conditions, when dead or 
crude oil is used, each pressure stage is prolonged for the time required to 
obtain wettability stabilization at a given saturation. Drainage is thus carried out 
in approximately two to three weeks, which typically corresponds to the kinetics 
of the chemical processes linked with aging. The simplicity of this method lies in 
the fact that it just requires real-time recording and plotting of the mean 
saturation and resistance to obtain a continuous resistivity index curve free of 
any artefact. The key point is due to the fact that the radial geometry of the 
electrodes allows examining the whole volume of the sample and to 
compensate for the non-uniform saturation profile which appears in the absence 
of capillary equilibrium. 

[0020] Under ambient conditions, the radial-electrode cell described for 
example in the aforementioned patent EP-A-974,839 is used. The complex 
impedance measurements were carried out at a fixed frequency of 1 kHz and 
the real part was extracted to calculate the resistivity index. The highest 
capillary pressure applied was 12 bars (for an interfacial tension y=35 mN/m ; 
under reservoir conditions, with crude oil, the maximum capillary pressure was 
reduced proportionally to the interfacial tension). 

Conductivity models 

Double-porosity conductivity model (DPC) 

[0021] The goal is to explain the bending of curves Rl at low saturation 
(type I, Fig. 1). The double-porosity model considered here is basically very 
close to those proposed for clayey sands, where the clays present at the 
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surface of the pores constitute a parallel path for the current. For carbonates, 
we assume the existence of two pore networks having parallel electrical 
conductivities is assumed. The two main ingredients in the model are the 
description of the pore network invasion during drainage and the description of 
the electrical layout of the two different pore populations. 

[0022] The first network 1 (macropores for example) represents the major 
part of the pore volume, whereas the second network 2 (micropores) only 
represents a small fraction thereof, which is not necessarily greater than the 
percolation thresholds. The saturation of each network, Swi and Sw 2 , which are 
related to the measured mean saturation Sw is first considered by : 

Sw = f^Wi + f 2 Sw 2 where fi + f 2 = 1 (1 ). 

[0023] The quantities fi and f 2 represent the pore volume fraction of each 
population. These fractions were evaluated by means of the NMR relaxometry 
technique. We then assume that the networks are invaded by oil at different 
capillary pressures ; the oil reaches the small pores at a higher pressure than 
the pressure observed for the larger pores. From the capillary pressure curve 
previously drawn, this pressure corresponds to a mean saturation Sc which can 
be deduced from the curves obtained by mercury injection (showing the pore 
throat distribution). Sw-j can be expressed as a function of Sw at high saturation 
Sw : 

S Wl = Jl — -,Sw 2 =l for Sw>Sc (2) 
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[0024] Below Sc, relations Sw-i=f(Sw) and Sw 2 =f(Sw) require other 
hypotheses. The assumption (i) is a linear relation and (ii) that Sna^i ->0, Sw 2 ->0 
when Sw -^0. 

[0025] It is deduced that : 



fiSc 

Sw 2 = — for Sw < Sc (3) 
Sc 

[0026] The conductivity of each network is now considered. In case of an 
initial brine saturation, the total conductivity Ct 0 for the two networks in parallel 
will be : 

Ct 0 = Ci + C 2 = Ci (1 + a) where C 2 = a Ci (4). 

[0027] Parameter a is the conductivity ratio of the two 100 % saturated 
networks. From Archie's first law, Ci and C 2 are assumed to be related to the 
pore volume fraction of each population and, therefore, that a is of the order of 
ratio f 2 /fi. When the two networks are invaded by oil, each conductivity is 
assumed to be related to the saturation by a power law (as in Archie's second 
law). The total conductivity in the two saturation ranges is as follows : 

Ct = SWi n1 Ci + C 2 for Sw > Sc (5) 

Ct = Sw^ 1 Ci + Sw 2 n2 C 2 for Sw < Sc (6). 
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[0028] By means of equations 4, 5 and 6, the resistivity index Rl is as 
follows : 




1 + a Sw. 



1 + a 



— for Sw > Sc 

n\ J 



(7) 



RI = Sw, 



1 + a 



— for Sw < Sc 



(8) 



i + a Sw; n] I Sw; 



[0029] When a pore population is in the dominant position (network 1 in the 
present case), functions Rl are governed by the resistivity properties (ni) of this 
population. Equations 7 and 8 are essentially similar to the formulas used for 
clayey sands, except that a second exponent n 2 is introduced which 
characterizes the second network. 

[0030] In the DPC model, there are in total 4 parameters n-j, n 2 , a and Sc. 
For a given experimental curve, Sc is measured separately during a mercury 
injection experiment, whereas the other parameters are adjusted. A variation 
range is however observed for a around f 2 /fi, for which a physical explanation 
exists. At high saturation Sw, the slope on the log-log scale of RI(Sw) is -ni 
and, at low saturation, the slope is -n 2 . To a certain extent, Sc and a 
compensate for one another (Fig. 3), but a is the most sensitive parameter 
which controls the final value of Rl. It can be noted that the case Sc=0.05 
presented in Fig. 3 corresponds to a situation where the second network is not 
invaded by oil, which gives a horizontal asymptote. 
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[0031] The parameters of the model were adjusted in order to correspond to 
the experimental curve Rl measured on sample RC. The value of Sc=0.4 was 
deduced from the elbow of the capillary pressure curve obtained by mercury 
injection and the volume fraction fi=0.88 macropores (network 1) was deduced 
by NMR (Fig. 2). Despite the low permeability of the sample, the lowest 
saturation obtained is very low (4 %), which allows good determination of the 
parameters of the model under ambient conditions. The initial slope of curve Rl 
is ni=1.71 ; it characterizes network 1. Network 2 is very weakly sensitive to 
saturation (n 2 =0.25) and does not behave like a standard network. However, the 
initial conductivity ratio <x=0.054 of network 1 to 2 is of the order of f 2 /fi=0.136. 
Under reservoir conditions (oil wettability conditions), a high increase of ^ is 
observed, but curve Rl is still not linear on the log-log scale. Characteristics 
n 2 =0.78 and a=0.11 of the second network are also slightly modified, but the 
precision concerning these parameters is lower that under ambient conditions 
because the saturation reached (at the same capillary pressure) is much higher 
and close to Sc. 

Approximate formula in the case of the DPC model 

[0032] An approximate formula reducing relations 7 and 8 to a single 
formula, which is valid for the type I curves is as follows. Note that : 

Sw 2 n2 Sw^ 1 « Sc n2 Sw n2 ' n1 » Sw n2 " n1 

with Swi = Sw and Sc n2 « 1 (15) 

which gives f 2 «f 1 . Since n 2 «ni, equations 7 and 8 can be approximated by : 

Rl = Sw~ nl 1 + C 2 : (16). 
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[0033] Using the 3 parameters rv» f ri2 and C is sufficient to describe the data 
with precision by means of equation 16. Here, the meaning of C is 
approximately the same as a and it can depend on the temperature. The 
formulation given in equation 16 can be used to describe standard experiments 
carried out at capillary equilibrium where only a limited number of points is 
available. 

Triple-porosity conductivity model (TPC) 

[0034] The goal is to explain the straightening of curves Rl at intermediate or 
high saturation, and their bending at low saturation on the log-log scale (types II 
and III, Fig.1). Complex carbonates can have three pore populations referred to, 
for simplification reasons, as micro, macro and mesopores (3, 2 and 1 
respectively). As in the double-porosity model, the saturation of the 3 
populations is considered: 

Sw = ^Swi + f 2 Sw 2 + f3Sw 3 where fi + f 2 + h = 1 (9) 

[0035] Here again, the invasion of these populations by oil during drainage is 
assumed to be sequential. If network 1 is invaded first, a mean saturation Sm is 
defined at which network 2 is invaded : 

s „ | = Wa-/> , 

Sw 2 =l Sw 3 =l for Sw > Sm (10) 
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[0036] Below Sm, it is possible to imagine many scenarios. In general, linear 
relations for functions Swi(Sw) and Sw2(Sw) are assumed. According to a 
possible scenario, Swi ->0[,] and Sw 2 ->0 when Sw-^Sc is assumed. Sc is the 
saturation at which the micropores are invaded. It follows therefrom that : 



[0037] This invasion scenario is summarized in Fig. 5. Typically, a situation 
is considered where f3«fi and fi«f2, and the micropores (network 3) are 
invaded at a pressure that is much too high to be observed during the 
experiment. 

[0038] Consider now the conductivity of these populations. Networks 1 and 2 
are assumed in series, whereas network 3 is in parallel as shown in Figure 6. 
The arrangement in series may seem to contradict a sequential invasion by oil 
(which is an essentially parallel mechanism), but such a situation is possible. At 
Sw=1 , conductivity Ct 0 of the system shown in Fig. 6 is as follows : 




Sw-Sc 
Sm-Sc 9 



Sw 3 =1 for Sc<Sw<Sm 



(11) 



Ct 0 = (CV 1 + CzY + C 3 = Ci [(1 + en" 1 )" 1 + a 2 ] 



where C2 = oci C1, C3 = 0.2 C1 



(12). 



[0039] In the two saturations ranges, it is obtained that : 



Ct = C1 [(Swr n1 + on" 1 )" 1 + a 2 ] for Sw > Sm 



(13) 



Ct = C1 [Sw^ 1 + Sw 2 " n2 af 1 )* 1 + a 2 ] for Sc < Sw < Sm (14). 
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[0040] The resistivity index can be calculated from equations 12 to 14. For 
the measured values of fi, f 2 and f3, n 1f n 2 , <xi and a 2 have to be adjusted to the 
experimental data. As in the DPC model, on is expected to be of the order of 
f 2 /fi and a 2 of the order of f^f 1 . 

[0041] This model was tested with resistivity index curves measured on 
samples EL and BL. On sample EL (Figure 7), straightening of the curve at high 
saturation is reproduced qualitatively at high saturation and flattening at low 
saturation is reproduced. A discontinuity at Sm (deduced from the elbow of the 
capillary pressure curve obtained by mercury injection) is observed because of 
a sudden modification of relation Swi(Sw) (as shown in Fig. 5). After 
adjustment, the two dominant pore populations have similar characteristics 
(ni=n 2 =1.5). The third population, which behaves like a parallel circuit, is not 
invaded by oil and it is only characterized by conductivity ratio oc 2 . On sample 
BL (Fig. 8), curve Rl can be explained by a similar mechanism of resistance in 
series, but with a very high value Sm. The two networks observed have very 
different values n and that network 2 is much less conducting than network 1 
(ai=7.1). 

[0042] The method according to the invention finds applications notably in 
the sphere of hydrocarbon reservoir development. 
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